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I.   INTRODUCTION 


A.   EACKGSCUKE 


As  cf  1S74,  the  United  States  contained  only  6%  cf  the 
world  population  but  was  using  33X  of  the  energy  consumed 
each  year  according  to  Bef.  1.  Approximately  one  third  of 
the  oil  used  in  the  U.S.  is  imported.  This  heavy  dependence 
en  foreign  sources  fcr  energy  coupled  with  the  exhaustion  of 
domestic  fcssile  fuels,  particularly  oil,  in  the  foreseeable 
future  is  causing  the  U.S.  military  and  civilian  sectors  of 
the  economy  tc  search  for  energy  sources  that  would  be 
inexhaustible  and  independent  of  foreign  control. 
Increasing  fuel  costs  and  environmental  problems  are  other 
forces  pushing  the  search  for  alternatives.  Nuclear  power, 
once  considered  the  cure-all,  is  facing  severe  problems. 
Costs  cf  construction,  operation,  and  fuel  have  risen 
dramatically.  Questions  concerning  the  safety  of  operation 
of  the  nuclear  plants  and  of  the  storage  cf  nuclear  wastes 
have  beccire  political  issues  that  have  slowed  the 
exploitation  cf  nuclear  power.  The  present  nuclear  reactors 
use  a  fuel  ttat,  like  petroleum,  has  a  limited  availatility 
and  whose  price  is  increasing  rapidly.  The  breeder  reactor, 
which  would  create  more  fuel  than  it  burns,  and  the  fusion 
reactor,  which  would  use  the  hydrogen  found  in  water,  are 
still  a  long  way  frcm  producing  power  for  a  world  whose 
demands  fcr  energy  increase  each  year. 

Many  pecple  advocate  the  exploitation  cf  the  so-called 
"free"   enercy  scurces.   These  sources  are:   (1)  the  kinetic 
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energy  of  moving  fluids,  such  as  winds  and  ccean  currents; 
(2)  the  potential  energy  of  tides  and  rivers;  (3)  the  heat 
generated  inside  the  earth;  (4)  the  direct  conversion  of 
solar  energy  into  electricity  and  heat;  and  (5)  the  use  of 
sun  warmed  water  in  conjunction  with  a  source  of  colder 
water  to  provide  the  temperature  differential  to  run  a  vapor 
power  cycle.  Bene  of  these  sources  are  truly  "free".  iihat 
is  happening  is  the  trading  off  of  the  transportation, 
processing,  and  environmental  costs  of  a  concentrated  energy 
source  for  the  capital,  operating,  and  social  costs  of 
converting  a  very  diffuse  source  into  a  mere  concentrated 
form  such  as  electricity. 

The  earth  can  te  thought  of  as  a  giant  heat  engine, 
absorbing  energy  from  the  sun  and  reflecting  energy  as 
thermal  radiation  tack  to  space  as  pictured  in  Hef.  2. 
Since  the  Ecles  receive  less  energy  than  the  equator,  the 
atmosphere  and  ocean  attempt  to  distribute  the  energy  more 
evenly  ever  the  earth.  The  air  is  heated  by  absorption  of 
sunlight  and  by  contact  with  the  surface  of  the  earth.  The 
heated  air  rises  thus  setting  up  surface  wind  currents  as 
cooler  air  tries  tc  replace  the  rising  air.  The  shear 
stress  created  by  the  relative  moticn  between  the  air  and 
the  surface  water  causes  the  water  to  move  thus  creating  the 
surface  water  currents.  Many  currents  in  the  oceans, 
regardless  cf  depth,  originate  with  the  wind  shear  stress  at 
the  air-surface  interface.  The  circulation  patterns  cf  the 
oceans  become  extremely  complex  due  to  the  influence  cf  the 
rotation  cf  the  earth  and  the  shapes  of  the  ccean  basins. 

In  general  the  surface  waters  near  the  equator  are 
warmed  by  the  sun  and  flew  towards  the  poles  giving  up 
energy  by  evaporation,  and  by  convection  tc  the  air  ard  by 
radiation  back  to  space.  As  these  currents  cool,  seme  of 
the  water  becomes  dense  enough  tc  sink  and  then  flews  back 
to  the  eguatcr  alcng  the  bottom. 
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The  relationship  between  salinity  and  temperature  and 
the  density  cf  seawater  is  such  as  tc  create  a  stable 
stratification  of  the  ocean  in  many  areas.  A  salinity 
decrease  or  a  temperature  increase  reduces  the  density.  The 
total  effect  is  to  create  a  warm  surface  current  overlaying 
the  cold  water  returning  from  the  poles.  In  the  tropical 
waters  the  surface  temperature  may  exceed  80°F  (26.7°C) 
while  3,Q0C  to  6,000  ft.  (914  to  1828  m)  dcwr  the 
temperature  nay  be  about  40°F  ^.4°C).  Ibis  temperature 
difference  can  be  used  to  rur  a  man-made  heat  engine,  like 
the  Ocean  Thermal  Energy  Converion  (OTJ.C)  power  plant 
ccrcept. 

The  pioper  choice  of  a  site  is  critical  to  the  cost  of 
such  a  heat  engine.  Since  the  temperature  difference  is  low 
ccapared  to  conventional  thermal  power  generating  methods 
which  typically  use  temperature  differentials  exceeding 
50C°F  (26C°C),  the  effect  of  a  1°F  (0.56°C)  loss  in  the  AT 
is  much  nore  costly  at  the  low  thermal  differentials,  which 
are  usually  less  than  50°F  (27.8°C)  in  the  ccean. 
Therefore,  the  vertical  temperature  profile  of  the  cceans 
will  ccntrcl  the  location  of  an  OTEC  power  plant.  The 
surface  temperatures  vary  with  location  and  with  time.  At 
mcst  places  in  the  ocean  there  is  a  seasonal  change  of  the 
surface  sater  temperature  of  5  to  8°F  (2.8  to  4.4°C)  except 
in  the  trcpical  areas  that  remain  relatively  constant 
according  to  Bef.  3.  However,  all  regions  cf  the  oceans  are 
affected  ty  wind-generated  mixing  of  the  surface  water  that 
changes  the  surface  temperature  by  a  degree  or  so  and  a 
daily  cycle  cf  a  few  tenths  cf  a  degree.  Some  cf  the 
possible  O.S.  near-shore  sites  being  considered  are  the 
waters  arcund  the  Eawaian  Islands  and  the  Gulf  Streai  off 
the  southeastern  U.S. 

In  the  past,  the  engineering  analysis  and  the  cost 
estimation   have   been   conducted   as   if  they  were  separate 
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functions.  In  a  system  with  a  small  temperature 
differential,  such  as  this  one,  the  engineering  design  and 
economic  analysis  should  be  linked  together  as  a  single 
procedure  because  of  the  capital  intensive  nature  cf  the 
ocean  thermal  power  plant.  For  example,  the  heat  transfer 
egiraticn  (C  =  DA  AT)  is  coupled  to  the  cost  equation  for 
heat  exchangers  ($  =  K  An)  through  the  surface  area  of  the 
tubes.  Eecause  cf  the  small  AT  involved  the  area  must  be 
very  large  tc  supply  the  same  Q.  The  other  components  are 
all  linked  together  by  the  engineering  and  cost  equations  so 
that  the  problem  cannot  logically  be  rroken  intc  an 
equipment  by  equipment  optimization.  The  present  analysis 
coosiders  a  closed  vapor  power  system  operating  cu  the 
thermal  potential  available  from  the  ocean.  The  economic 
optimization  cf  the  system  is  carried  out  to  shew  what 
information  can  he  gotten,  what  conclusions  can  be  reached 
frcm  such  analysis  and  where  research  effort  should  be 
expended  tc  improve  the  design  of  the  system.  This  study 
does  not  attempt  tc  predict  the  cost  of  such  a  system  cr  to 
make  specific  recoamendations  about  the  parameters. 

A  schematic  diagram  of  the  system  is  shewn  in  figure  1. 
The  warm  surface  water,  the  Gulf  Stream  for  example,  is 
pumped  through  a  heat  exchanger.  There  seme  energy  is 
transferred  ty  boiling  the  working  fluid  such  as  ammonia. 
The  working  fluid  is  piped  to  a  turbine  where  some  of  the 
energy  is  ccnyerted  to  mechanical  energy  then  to  electrical 
energy  ic  a  generator.  The  unavailable  portion  of  the 
absorbed  energy  is  rejected  to  the  environment  in  the 
condenser  and  the  working  fluid  is  returned  to  the  liquid 
state.  The  lew  temperature  in  the  condenser  is  maintained  by 
pumping  ccld  ocean  water  up  from  the  depths.  Finally  the 
working  fluid  is  pumped  back  to  the  boiler. 

The  Ccear  Thermal  Energy  Conversion  program  has  received 
much    attention    in   the   last   several   years   from   many 
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FIGURE  1.  SYSTEM  SCHEMATIC  DIAGRAM 
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researchers  funded  at  first  primarily  by  the  National 
Science  Fcurdation  (NSF) .  Since  the  establishment  cf  the 
Energy  Research  and  Development  Administration  (ERDA)  ,  new 
research  ccrtracts  are  funded  mainly  by  this  agency.  The 
C.S.  Navy  is  alsc  funding  some  research  in  the  OTEC  area. 
Bef.  4,  5  and  6  give  a  good  summary  of  the  work  that  has 
been  done.  ihile  much  of  the  work  has  concentrated  on  the 
heat  exchangers,  the  other  components  and  the  system  as  a 
whcle"  are  receiving  attention.  Some  researchers  are 
concentrating  on  the  environmental  and  legal  questions  of 
such  systems.  Cther  teams  are  considering  hew  to  use  the 
energy  supplied  from  an  OTEC  power  plant. 

In  any  design,  trade-offs  must  be  made.  Capital  cost 
can  be  traded  for  operational  costs.  Pumping  power  can  be 
traded  for  heat  exchanger  size.  Trade-offs  can  involve 
items  that  have  a  common  measure  such  as  dollars;  seme  other 
trade-offs  are  very  hard  to  measure  in  common  terms  such  as 
social  and  environmental  costs.  The  trade-effs  may  involve 
crly  one  cempenent  such  as  a  pump  or  a  turbine,  but  most 
often  the  trade-effs  affect  more  than  one  component  or 
possibly  the  entire  system. 

The  trade-off  process  should  continue  until  an  optimum 
design  is  obtained.  From  Ref.  7,  "The  aim  cf  optimization 
is  the  selection,  cut  of  the  multiplicity  of  potential 
solutions,  of  that  solution  which  is  the  best  with  respect 
to  some  well-defined  criterion."  In  order  for  the  plant 
design  tc  be  optimum,  the  components  must  be  optimized 
within  the  context  of  the  entire  system.  Some  cf  the 
parameters  for  a  piece  of  equipment  may  have  no  direct 
effect  en  the  optimum  design  of  another  portion  cf  the 
system,  but  seme  parameters  will  have  significant  effects. 
For  example,  a  heat  exchanger  is  made  up  cf  a  large  number 
of  tubes  cf  a  given  diameter  and  length  and  is  designed  to 
have   a   certain  heat  transfer  rate  for  a  given  flow  rate  of 
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water  pumped  through  the  tubes.  The  diameters,  length  and 
the  number  cf  tubes  may  be  optimized  for  lowest  cost  cf  the 
heat  exchanger.  This  approach  leaves  out  the  effect  that 
the  tube  diameter  and  the  tube  length  have  en  the  pressure 
drop  across  -the  heat  exchanger.  Minimizing  the  cost  cf  the 
heat  exchanger  involves  minimizing  the  surface  area.  This 
may  lead  tc  increasing  the  pressure  drcp  across  the 
exchanger.  The  higher  the  pressure  drop  the  more  powerful 
and  the  mere  expensive  the  pump  must  be  to  maintain  the  flow 
velocity.  This  is  just  one  example  cf  the  hundreds  of 
trade-offs  possible. 

Unfortunately,  what  happens  in  many  cases  is  that  each 
group  of  experts  working  on  a  particular  component  optimize 
that  item  around  some  given  conditions  that  may  have  keen 
set  by  teams  working  en  some  other  section  cf  the  sjstem. 
Little  analysis  is  made  of  how  the  entire  system  respcEds  to 
the  coupling  constraints.  All  that  the  pump  manufacturer 
wants  is  that  the  buyer  tell  him  what  the  flow  rate,  head 
and  service  conditions  are  and  he  will  quote  a  price. 
Eefcre  the  buyer  can  decide  on  the  flow  rate  and  head  he 
should  first  knon  how  the  cost  cf  the  pump  varies  as  flow 
rate  and  head  vary  so  that  those  two  parameters  can  be 
optimized  in  the  context  of  the  system. 

The  prcclem  then  becomes  one  of  acquiring  sufficient 
information  en  hew  costs  vary  as  certain  parameters  change. 
Each  manufacturer  is  able  to  make  estimates  cf  costs  for  his 
product  but  he  is  reluctant  to  devulge  his  information  to 
buyers  or  cempetitcrs.  This  makes  constructing  cost  curves 
a  hazardous  prccess  at  best,  since  the  data  must  be  gathered 
by  ether  nettcds  that  may  cause  large  unexplained  variances. 

Sometimes  the  optimization  analysis  is  made  with 
assumptiens  that  make  the  answers  of  questionable  value. 
For  example,  sometimes  the  heat  transfer  coefficient,  D,   is 
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assumed  fixed  in  the  equation  g  =  0A  AT.  Eut,  as  shewn  in 
Ref.  8,  for  boiling  the  heat  transfer  coefficient  is  very 
sensitive  to  the  temperature  difference  with  the  result  that 
g  ■  UA,A13«33.  lhe  The  accuracy  of  the  engineering  eguations 
should  re  kept  in  mind.  Most  of  the  equations  are 
ccrrelaticns  tc  fit  experimantal  data.  The  correlation  was 
developed  fcr  seme  particular  set  of  data  and  may  disagree 
significantly  when  compared  to  data  taken  by  other 
experimenters.  The  limitations  of  the  various  relations 
must  be  kept  in  mind  when  reviewing  the  results.  For  a 
system  the  size  of  the  OTEC  plant,  a  pilot  plant  should  be 
constructed  in  order  to  prove  out  the  answers  given  by  the 
calculations  before  expensive  mistakes  are  made. 

As  a  design  progresses,  the  analysis  begins  with  a 
general  picture  containing  many  simplifying  assumptions  and 
prcceeds  foreward  with  more  details  added  at  each  stage. 
This  report  is  the  second  one  written  at  the  Naval 
Pcstraduate  £chccl  on  the  subject  of  OTEC.  Commander  Furman 
Sheppard,  OSH,  in  his  analysis,  Ref.  9,  attempted  tc  show 
what  kinds  cf  information  could  te  gained  from  a  combined 
thermal  econcmic  model.  He  advocated  the  use  cf  zone 
analysis,  as  developed  by  Wismer  in  Ref.  10,  fcr  large 
complex  systems.  In  Wismer* s  approach  to  zone  analysis,  the 
system  is  brcken  down  into  zones  containing  one  cr  more 
components .  Ihe  zones  are  connected  by  linking  variables. 
Two  methods  cf  optimization  are  possible.  Fcr  the  first 
method,  in  each  zone,  the  zone  parameters  are  optimized 
while  the  Hiking  variables  are  held  fixed.  Then,  the  zone 
parameters  are  held  fixed  while  the  linking  variables  are 
optimized.  lie  process  iterates  between  optimizing  the  zone 
parameters  and  the  linking  variables  until  convergence 
criteria  are  satisfied.  In  the  second  method  the  zenes  are 
cut  off  frcii  each  other  and  each  linking  variable  is 
separated  intc  two  variables,  one  on  eacb  side  of  the 
cutting   flare.   The  optimization  proceeds  until  the  linking 
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variables  en  each  side  of  the  cuts  are  equal. 

Sheppard»s  system  did  not  produce  vapor  in  a  boiler  but, 
instead,  maintained  the  working  fluid  in  the  liguid  state  in 
the  warn  seawater  heat  exchanger  and  vaporized  the  working 
fluid,  aittcnia  in  his  model,  by  expanding  it  through  a 
"black  fccx"  called  a  vaporizer.  His  optimization  analysis 
was  confined  to  the  zone  containing  the  warm  seawater 
exchanger,  the  feed  pump,  and  the  warm  water  circulating 
pump. 


E.   CEJECIIVIS 


The  present  analysis  improves  the  previous  model 
analysis  ty  replacing  the  liguid-to-liguid  heat  exchanger 
and  the  vaporizer  with  a  boiler.  The  boiler  model  is  to  be 
of  sufficiert  detail  so  that  a  realistic  design  is 
represented.  The  ether  uajor  components  of  the  system,  such 
as  the  turbice,  the  condenser,  and  the  condenser  circulating 
pumps,  are  included.  A  nonlinear  programming  technigue  is 
used  to  perform  the  thermal-economic  optimization  of  the 
completed  system.  The  research  reported  here  is  intended  to 
show  the  csefulness  of  optimization  analysis  in  the  design 
of  a  lew  thermal  difference  system  and  to  demonstrate  some 
cf  the  prcblens  that  must  be  considered  in  this  type  of 
analysis. 
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II.   THE  MODEL 


A.   GENES2L  EFSCBIfTION 


Vapcr  power  generating  cycles  usually  consist  of  four 
processes.  Energy  is  absorbed  from  a  high  temperature 
source  and  is  used  to  vaporize  a  working  fluid  in  a  bciler. 
Next,  work  is  extracted  by  expanding  the  vapcr  in  an  engine, 
exhaustinc  the  vapcr  at  a  lower  temperature  and  pressure. 
The  working  fluid  is  then  condensed  back  to  the  liguid  state 
by  rejectirg  energy  to  the  lower  temperature  energy  sink  in 
the  third  prccess.  Finally,  the  liquid  is  returned  to  the 
higher  pressure  of  the  boiler  by  a  pump.  In  the  following 
analysis,  the  energy  source  is  a  current  of  ocean  water 
heated  bj  the  sun  and  the  energy  sink  is  assumed,  to  be 
colder  ocean  water  from  deep  currents.  The  ticrking  fluid  is 
assumed  to  be  propane  although  many  other  fluids  are 
possible,  such  as  ammonia  or  the  Ereons. 

Six  ccmpcnents  are  modeled  in  the  following  analysis. 
The  hot  side  heat  exchangers,  the  boilers,  vapcri2e  the 
working  fluid  which  is  expanded  through  a  turbine  exhausting 
to  a  condenser  that  rejects  the  unused  heat  to  the  cold 
sink.  Three  pumps  are  included,  a  feed  puirp  returns  the 
working  fluid  to  the  boiler  and  two  pumps  circulate  the  warm 
and  cold  waters  through  the  bciler  and  the  condenser, 
respectively.  The  arrangement  of  the  system  is  shown  in  Fig. 
2.  The  nodes  are  numbered  to  correspond  with  the  fluid 
state  points  shown  later,  in  Fig.  3. 
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The  analytical  model  consists  of  the  equations 
describing  the  system  and  a  set  cf  parameters  to  be 
optimized  sc  as  tc  minimize  the  capital  cost  of  the  plant 
subject  tc  certain  constraints.  The  objective  function,  the 
function  tc  te  minimized,  is  the  sum  of  the  capital  costs  of 
the  pumps,  the  boiler  and  the  condenser.  Although  the 
turbine  is  a  significant  capital  cost  and  would  affect  the 
design  cf  the  system,  it  is  not  included  in  this  model. 
There  are  twc  classes  of  constraints.  The  first  class  is  the 
explicit  constraints  that  specify  the  upper  and  lower  bounds 
on  the  parameters.  The  other  class  of  constraints  is  the 
implicit  constraints.  These  are  the  engineering  eguations 
that  describe  the  operation  of  the  power  plant.  For 
example,  a  constraint  might  say  that  the  temperature  inside 
the  boiler  mcst  be  less  than  the  temperature  of  the  warm 
seawater  because  the  eguations  fail  if  the  requirement  is 
net  satisfied.  Another  example  would  be  a  power  balance 
that  reguires  that  the  sum  of  the  energy  flews  into  and  out 
of  the  system  be  zero.  Additionally,  they  proscrib€  the 
region  that  the  optimization  rcutine  may  search  within  for 
the  minimum  cf  the  problem.  The  implicit  constraints  are  of 
twc  types,  eguality  and  inequality  constraints.  Ineguality 
constraints  set  one-sided  bounds  on  the  feasible  region. 
Seme  of  the  ineguality  constraints  are  restrictions  placed 
on  the  maximum  or  minimum  values  of  the  parameters  dee  to 
judgements  Bade  en  the  basis  of  practical  engineering 
consideraticrs.  The  equality  constraints  are  relationships 
that  must  be  strictly  satisfied  at  the  optimimum.  The  system 
of  eguations  describing  the  model  fit  into  the  following 
form: 
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nirioiize  f  (X) 


subject   to:      g    (X)    >   0,    i=1,2,....,m 
i 


h    (X)    =    0,    i=m+1 , . . . . ,m+n 
i 


where  X  is  the  vector  of  the  parameters  tc  be  optimized, 
f(X)  is  the  capital  cost  function,  g^  (X)  are  the  inequality 
constraints  and  hj  (X)  are  the  eguality  constraints.  The 
parameters  are  the  dimensions  of  the  toilers  and  the 
condensers,  the  boiling  and  condensing  pressures,  the  sea 
water  velocities,  and  the  mass  flew  rate  of  the  wcrking 
fluid  through  the  system. 

The  tasis  for  the  economic  framework  must  be  set  cut  at 
the  start.  All  measurable  costs  are  relevant  to  the  final 
design  cf  the  power  plant,  but,  for  some  types  of  analysis 
and  design,  iiany  ccsts  can  be  ignored.  The  life  cycle  cost 
of  the  C1EC  must  be  considered  before  the  decision  is  made 
to  commit  the  D.S.  and  the  rest  of  the  world  tc  the 
censeguences  cf  choosing  this  energy  generation  scheme  as 
even  a  partial  solution  to  the  energy  problem.  The  life 
cycle  ccst  consists  of  the  capital  costs,  the  financing 
costs,  and  the  operating  costs.  The  trade-off  of  operating 
ccst  and  capital  cost  is  not  considered  in  this  research, 
because  the  operating  cost  is,  for  the  most  part, 
independent  cf  the  overall,  dimensional  characteristics  of 
the  components.  Operating  costs  are  dominated  by  the 
structure,  site  location,  and  maintenance  requireuents. 
Such  things  as  the  kind  of  instrumentation,  controls, 
bearings,  and  structural  materials  have  mere  effect  en  the 
operational  ccsts  than  the  length  of  the  bciler  or  the 
diameter  cf  the  toiler  shell.  Therefore,  enly  the  capital 
ccsts  of  the  majcr  components,  except  the  turbine,  are 
considered. 
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To  b€  cseful  for  model  analysis  cost  data  must  be 
transformed  into  a  cost  estimating  relationship  (CEB)  that 
shews  hew  cost  varies  with  some  variable  or  some  combination 
of  variables.  The  costs  must  be  demonstrated  to  be  highly 
correlated  with  the  variable  chosen.  Since  the  cost 
relation  is  enly  a  correlation,  the  bounds  on  the  accuracy 
of  the  eguaticc  should  be  specified  by  those  who  develop  it. 
Since  the  relation  is  arrived  at  from  a  finite  body  cf  data, 
the  CEB  is  directly  applicable  only  in  the  range  of  the  data 
collected.  Extrapolation  cf  the  CEB  outside  the  range  is 
dangerous  because  there  is  no  data  in  that  region  to  support 
the  assumed  curve.  However,  the  value  of  CEB*s  is  the  use 
of  them  to  enable  the  estimator  to  make  a  better  analysis  of 
a  new  design  than  he  could  if  the  information  is  not 
available.  The  relevant  range  of  the  CEB  must  be  kept  in 
■ind  when  analyzing  the  solutions  from  the  standpoint  cf  the 
reliability  cf  the  answers.  Another  factor  to  consider  when 
using  cost  estimating  relationships  is  whether  the  equipment 
the  data  was  taken  from  is  comparable  to  the  systen  the 
analyst  is  considering.  If  the  design,  the  technology,  or 
the  application  of  the  proposed  equipment  is  significantly 
different  ficm  that  from  which  the  data  is  taken,  the  CEB 
may  give  unreliable  answers. 

This  analysis  begins  with  the  development  of  the 
engineering  eguatiens  that  constrain  the  system.  This  is 
followed  by  the  development  of  a  cost  framework  based  on 
capital  costs  as  the  function  to  be  minimized.  The 
engineering  and  economic  relations  are  linked  together  in 
the  objective  function  and  the  constraints.  A  description 
of  the  algorithm  that  performs  the  optimization  is  included. 


B.   CYCLE  AN2IYSIS 
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The  building  of  the  model  begins  with  the  analysis  of 
th€  therncdynamic  cycle.  Figure  3  is  a  temperature-entropy 
diagram  of  the  working  fluid  on  which  the  cycle  is  shown. 
These  numbered  nodes  represent  the  thermodynamic  state 
points  at  that  location  in  the  system  and  correspond  tc  the 
numbers  ic  fig.  2.  The  boiling  and  condensing  processes  are 
assumed  tc  be  reversible  processes.  In  the  pump  and  the 
turbine  the  reversible  process  is  indicated  by  the  subscript 
"s"  on  the  state  pcint  number.  In  the  actual  system  there 
are  irreversible  losses  inside  each  component  that  are  left 
out  of  the  analytical  model  or  lumped  in  with  ether  lesses. 
The  pressure  drops  due  to  frictional  losses  in  the  toiler 
and  the  cendenser  are  neglected. 

State  pcint  1  assumes  that  the  liquid  leaving  the  toiler 
is  saturated.  State  points  2s  and  2  represent  the 
conditions  leaving  the  feed  pump.  If  there  were  no  lesses  in 
the  pump,  the  fluid  would  be  at  2s.  "  The  fluid  enters  the 
bciler  tc  be  heated,  first,  to  saturation  conditiens  at 
state  point  2  then  boiled  at  essentially  constant  pressure. 
State  pcint  4  shows  the  vapor  leaving  the  bciler  as  a 
saturated  vapor.  The  lines  from  point  4  to  5s  and  5 
represent  tte  expansion  through  the  turbine.  The  path  tc 
state  pcint  5  accounts  for  the  losses  due  tc  frictior  and 
ether  irreversible  processes  in  the  turbine.  The  remaining 
heat  is  then  rejected  from  pcint  5  to  point  1,  again  assumed 
to  be  a  censtant  pressure  process. 

The  ttecretical  beat  transfer  rate  in  the  boiler  and  the 
condenser  is  found  by  subtracting  the  enthalpy  at  state 
points  2  and  1  from  the  enthalpy  of  state  points  4  acd  5, 
respectively,  and  multiplying  by  the  total  mass  flow  rate  of 
the  working  fluid. 

Qab  ■  M^(fl4  -  Ha)  (1) 
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FIGURE  3.  THERMODYNAMIC  CYCLE 
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The  enthalpy  at  state  point  5  is  determined  by  first  finding 
the  moisture  content  of  the  fluid  exiting  the  turbine. 


i  --2-i.  (3, 

s'9 


The  enthalpy  at  point  5s  is  determined  from 


HSs  =  H1  +  i  H,g  <4) 


This  aliens  E5  to  te  found  from 


?tur, -!aj_!s*  ,S) 

where   the   efficiency  of  the  turbine  is  assumed  kncwc  The 
power  out  of  the  turbine  is  nou  known. 


ptur  ■  *W<HM  -  HS)  <6) 

For  the  feed  pump  the  change  in  enthalpy  is  accounted  for  by 
tbe  change  ir  pressure  assuming  negligible  rise  is  the 
temperature.  The  enthalpy  of  point  2  is  determined  frca 

?     .  '1(P»5  "Pi)  (7) 

P      H2  "  H1  . 

where  the  efficiency  is  assumed  known.  The  power  reguired  to 
drive  the  juop  is  determined  from 
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*ip     =    Mw{<H2    "    H1>  <8> 


Tc  gain  a  better  feeling  fcr  the  low  thermal  potential 
systems  and  tc  make  comparisons  with  other  systems  several 
efficiencies  are  worthwhile  calculating.  The  first  is  the 
Carnot  efficiency  which  defines  the  maximum  efficiency 
possible  for  a  system  operating  between  two  temperature 
liiits  ic  decrees  absolute/  °H  (°K). 


\»r    =  1  >  ^2i±.  <9) 

Thot 


The  next  efficiency  calculation  takes  into  account  the 
theoretical  thermodynamic  cycle  and  is  the  ratio<of  the  net 
poser  out  of  the  cycle  to  the  heat  taken  up  by  the  cycle 


(eye  - 


_    gtur  *  *ip 


(10) 


This  acccimts  for  the  thermodynamic  inefficiencies  inherent 
in  the  cycle.  Finally,  the  system  efficiency  takes  into 
acccunt  the  ether  internal  consumers  of  power:  the 
circulation  pumps,  the  pipe  losses,  and  any  other  lesses, 
such  as  the  antifouling  eguipment.  This  calculation 
accounts  fcr  all  of  the  mechanical  and  electrical 
ineff iciercies  of  the  system. 


„      NET  POWEB  AT  THE  EUS  EfiB 

?sys    =  ■ (11) 

Qab> 


For    example,    if    the   hot    water   is    at      75°F       (23.9°C)      and 
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the  ccld  water  is  assumed  at  40°F  (4.4°C)  the  Carnot 
efficiency  is  6.5ft.  Fcr  one  typical  set  of  assumptions,  as 
shewn  in  Table  5,  where  the  turbine  and  the  pumps  are 
assumed  tc  be  85X  efficient,  the  cycle  efficiency  is  2.69? 
and   the  system  efficiency  is  2.455   for  one  feasible  design. 


C.   ENGIKEEBIHG  FRAMEWORK 


This  section  describes  the  design  of  the  components  that 
make  up  the  system.  The  level  of  detail  of  the  analysis 
provides  a  physically  realistic  model  and  clcses  the  lcop. 
Tte  pipe  friction  losses  on  the  cold  and  hct  water  inlets  as 
well  as  tte  pipe  losses  in  the  working  fluid  side  are 
neglected.  Ihe  turbine  is  included  only  as  an  energy 
transducer  and  no  attempt  is  made  to  obtain  costs,  physical 
dimensions,  cr  fluid  flow  characteristics.  The  efficiencies 
of  the  turbine  and  the  pumps  are  assumed  as  known  constants. 
The  heat  arscrbed  in  heating  the  working  fluid  to  the 
saturated  state  is  a  small  portion  of  the  total  energy  taken 
up  in  the  boiler;  therefore,  the  sensible  heating  of  the 
fluid  is  neglected. 

1 .   The  Eciler 

The  bciler  consists  of  a  number  of  tube-in-shell 
heat  exchangers  connected  in  parallel  as  shewn  in  Fie.  2. 
Figure  4  shoss  the  internal  arrangement  of  one  shell.  The 
tube  bundle  is  horizontal  and  boiling  occurs  on  the  outside 
surface  of  the  tubes.  The  tube  bundle  is  submerged  in  the 
working  fluic.  The  working  fluid  is  maintained  at  a  le^el  to 
just  cover  the  top  row  of  tubes.  The  warm  seawater  is  jumped 
through  the  inside  of  the  tubes. 
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FIGURE  4.  BOILER  CROSS  SECTION 
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The  clearance  between  the  tubes  is  fixed  at  0.5 
inches  (1.27  cm)  to  allow  for  flaring  of  the  tube  ends  or 
ether  space  and  strength  requirements.  Similarly,  the  tube 
thickness  is  fixed  at  0.028  inches  (0.071  cm),  since  no 
strength  requirements  due  to  differences  between  internal 
and  external  pressures  on  the  tubes  or  material  loss  due  to 
corrosion  of  the  tube  surfaces  are  accounted  for. 

At  the  tcp  of  the  tube  bank  the  leaving  vapor 
contains  a  large  amount  of  liquid  droplets  which  must  be 
eliminated  either  in  the  boiler  or  by  an  external  separator. 
In  the  mcdel  the  separation  takes  place  in  the  boiler  by 
allowing  sufficient  space  above  the  tube  bank.  Ealen  and 
Small,  Bef.  11,  presented  a  methed  for  finding  the  necessary 
separation  volume  that  considers  the  effects  of  surface 
tension,  density  of  the  vapor  and  the  liquid,  and  the  mass 
flew  of  the  jerking  fluid.  First  the  allowable  vapor  lead  is 
found  frcn 


r     o-       -i0-5 

VL  =  2290A  - (12) 

ri   l_6.  86x10-5  (^  -  p3)j 

Vapcr  load  is  an  expression  for  the  mass  flow  rate  cf  vapor 
through  a  ucit  volume.  The  expression  for  the  maximum  vapor 
lead  limits  the  flew  rate  so  that  the  moisture  droplets  are 
able  to  settle  out  of  the  vapor  before  it  leaves  the  bciler. 
TbeE  the  reguired  segmental  area  is  found  from 


SA  =  — (13) 

VL  L   S 


The  segmental  area  is  the  cross  sectional  area  from  the  top 
of  the  tube  bank  to  the  top  of  the  shell.  Kncwing  the 
segmental  area  the  height  of  the  tube  bank   above   the   main 
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diameter  is  found  from 


2       (2   hs\  2      2  0.5 

SA  =  C.25  Es   ARCOS J-  hs(0.25  Is      -  bs  )       (14) 


m- 


Palen  and  Small  also  present  a  much  simpler  method  of 
determining  the  total  height  of  the  tube  bank.  The  height 
of  the  separation  volume  is  set  at  4055  cf  shell  diameter 
which  results  in  "hs"  being  equal  to  0.1  Ds.  Either  aethod 
can  be  used  alone  or  the  first  method  can  be  used  as  a 
constraint  on  the  second.  For  this  problem  the  segmental 
area  is  fcund  by  bcth  methods  with  the  first  used  as  a  check 
on  the  values  given  ty  the  second. 

Since  the  heat  flux  calculation  is  based  on  a  single 
tube,  the  nunber  of  tubes  per  shell  must  be  determined  in 
order  tc  find  the  total  energy  absorbed  by  the  heat 
exchanger.  The  tubes  are  assumed  to  be  laid  out  in  a  60 
degree  triangular  pattern  for  maximum  compactness.  The 
tubes  are  located  at  the  vertices  of  an  eguilateral 
triangle.  Iwc  cross  sectional  areas  are  calculated.  The 
area  of  a  hexagon  whose  minor  radius  is  egual  to  one  half 
the  pitch  is  determined  from 
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2 

TAN  30°  (15) 


Then  the  effective  cross  sectional  area  of  the  tube  bank  is 
fcund  by  subtracting  the  clearance  area,  SA,  from  the  cross 
sectional  area  of  the  shell  after  correcting  the  shell 
diameter  for  a  clearance  space  to  allow  for  free  flow 
channels  both  within  the  tube  bundle  and  arcund  the  outside 
of  the  tuhe  bundle.  The  effective  cross  sectional  area  of 
the   tube   bark  is  divided  by  the  area  of  the  hexagon,  found 
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atcve,  tc  determine  the  number  of  tubes  in  a  shell. 

The  flew  of  the  liquid  and  vapor  through  the  tube 
bank  is  assuned  to  be  turbulent.  The  incoming  subcccled 
working  fluid  from  the  feed  pump  nixes  with  the 
recirculating  fluid  returning  frcm  higher  up  in  the  tube 
bank.  The  circulation  of  the  fluid  is  induced  ty  the 
boiling  process.  This  mixture  enters  the  tube  tank  and 
begins  the  toiling  process.  Near  the  bcttom  of  the  tube 
bark  the  teat  transfer  goes  mostly  to  bring  the  fluid  up  to 
the  boiling  temperature.  After  the  fluid  reaches  saturated 
liguid  conditions,  stable  boiling  begins.  The  vapor  bubbles 
rising  through  the  surrounding  liquid  increase  the 
turbulence  around  the  tubes.  As  the  liguid  continues  up 
through  the  tube  bank,  the  local  boiling  temperature  is 
decreasing  due  to  the  decreasing  hydrostatic  head.  This 
decrease  in  the  boiling  temperature  increases  the 
tenperature  difference  and  increases  the  heat  fie*  and 
increases  the  bubble  generation  rate.  Tne  vapor  generation 
may  become  sc  rapid  that  a  condition  called  vapor  blanketing 
may  ccccr  in  the  interior  of  the  tube  bank  where  the  liguid 
cannot  flow  inward  fast  enough  to  displace  the  vapor  flowing 
through  the  region.  This  is  similar  to  burnout  in 
cenvective  tciling  in  that  the  tubes  dry  out  and  the  heat 
transfer  rate  decreases.  This  effect  becomes  pronounced 
enly  near  the  critical  heat  flux  according  to  Starczewski, 
in  Hef.  8. 

All  cf  the  heat  and  nass  flow  variations  in  the  tube 
bundle  are  averaged  by  performing  the  calculations  for  a 
siEgle  tube  located  at  the  mid  height  of  the  tube  bank. 
This  tube  is  assumed  to  represent  the  average  of  all  the 
tubes  in  the  tube  tank.  The  boiling  situation  is  assuned  to 
be  saturated  pool  boiling.  Quoting  frcm  Collier,  Ref  12, 
"pool  boiling  is  defined  as  boiling  frcm  a  heated  surface 
submerged   in  a  large  volume  of  stagnant  liguid".   Since  the 
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fluid  volume  seen  by  the  tube  is  neither  large  nor  stagnant, 
a  correction  is  needed  for  the  effect  of  being  in  a  tube 
bundle.  Saturated  boiling  is  defined  as  boiling  at  a 
constant  temperature  and  pressure  as  the  liquid  changes  to  a 
vapor.  This  is  probably  not  strictly  true  because  the 
boiling  temperature  is  controlled  by  the  local  pressure  at 
€ach  tube,  and  the  fluid  is  not  reaching  a  stable 
temperature  as  it  rises  in  the  bundle.  The  situation  is 
probably  closer  to  convective  boiling  since  the  fluid  is 
boiling  in  a  confined  channel  between  the  tubes.  An 
additional  effect  in  the  boiler  is  the  fact  that  the 
hydrostatic  head  at  the  bottom  of  the  tube  bundle  causes  the 
boiling  tc  take  place  at  a  higher  pressure  at  the  tottom 
than  the  ceiling  at  the  top  of  the  bundle  where  the  pressure 
is  egual  tc  the  toiler  exit  pressure.  This  effect  is 
accounted  for  by  using  the  average  head  of  the  bundle  in  the 
heat  transfer  equations. 

The  heat  transfer  analysis  of  the  boiler  begins  with 
the  basic  equation 

Q  =  UA  ATm  (16) 

where  0  is  the  overall  heat  transfer  coefficient,  A  is  the 
area  through  which  the  heat  is  transferred  and  ATM  is  a 
prcperly  defined  mean  temperature  difference.  In  this 
paper,  the  product  UA  is  treated  as  a  single  guantity  called 
the  thermal  conductance.  One  may  look  at  the  equation  as  <m 
analogue  of  the  basis  electrical  eguation  I  =  E/R  where  I  is 
equivalent  tc  Q,  E  is  the  same  as  ATm  and  UA  is  1/E.  This 
analoqy  is  used  when  developing  the  expression  for  UA. 

In  order  to  account  for  the  effects  of  the  tube 
seeing  many  nearby  neighbors,  the  equation  is  modified  by 
the   inclusion  of  a  factor,  F,  to  account  for  the  effects  of 
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the  surrcundirg  tubes.   Equation  16  is  changed  to 

Q  =  OA  F  ATW  (17) 

where  F  is  defined  as  the  ratio  of  the  heat  flux  of  the  tube 
hurdle  tc  the  heat  flux  for  a  single  tube  in  a  stagnant  pool 
boiling  situation.  According  to  Ref.  13,  a  report  published 
by  workers  at  Heat  Transfer  Research,  Inc.  (HTRI)  of 
Alhambra,  California,  F  is  greater  than  one.  In  Figure  5, 
which  is  sketched  from  a  figure  in  their  report,  F  ranges 
frcn  4. 73  at  ATm  =  8°F  to  2.63  at  ATm  =  40°F.  Since  nc  data 
is  given  en  the  conditions  of  the  test,  F  is  assumed  to  be 
fixed  at  three. 

Ie  crder  to  find  AT,*,,  the  heat  exchanger 
effectiveness,  6  ,  is  determined.  The  effectiveness  is 
defined  as  the  ratio  of  the  actual  heat  transfer  rate  to  the 
maximum  ccssible  heat  transfer  rate  as  shown  by 


(18) 


where 


c  =  <M  c»>Hu,d  <19> 

C^m  =  minimum  of  C^  or  Cc        (20) 

In  the  bciler,  Th  m  is  the  temperature  of  the  entering 
seawater,  and  Tc In  is  the  saturation  tenperature  cf  the 
working  fluid  at  the  pressure  at  the  mid  height  of  the 
bciler.  Reference  14,  by  Kays  and  Lcndon,  contains 
eguations  for  finding  the  effectiveness  fcr  various  flow 
conf iguraticcs.    In   addition  to  finding  Cm;^ ,  CmQX  must  be 


36 


MAXIMUM                                             / 

HEAT   FLUX                                       / 

/~\ 

CM 

' 

+-> 

M-4 

11 

JZ 

P 

H 

» 

v_^> 

O* 

/              / 

o> 

/            47 

o 

/             0/ 

-4 

/            *7 

3/ 

^7 

^ 

O/ 

57 

&/ 

f 

*•*/ 

*V 

7 

f 

LOG  AT  (°F) 
FIGURE  5.  COMPARISON  OF  A  SINGLE  TUBE 
VERSUS  A  TUBE  IN  A  BUNDLE 


37 


kncwn  in  crder  to  decide  which  equation  applies.  For  the 
case  of  ceiling,  Cmax  is  infinite;  therefore,  regardless  of 
the  flew  pattern 


£  =  1  -  EXP(-NIO) 


(21) 


where   NTO,   the   number  of   heat  transfer  units  cf  a  heat 
exchanger,  is  defined  as 


NTD  =  UA/C^jn 


(22) 


and 


Cmin  ~  " 


iw    ^P.s>w 


(23) 


AT-is  found  frcm 


£  = 


OA  AT** 


(24) 


and    the    temperature   of   the   sea   water   at   the      outlet      cf      the 
tute    is    deternined   frcm 


£   = 


_    ThJ*    "    Th,out 


AT, 


(25) 


wfcere 


ATo    ~    Th,in~    Tc,in 


(26) 


Ibe   guantity   OA   is   determined   next.    Figure   6 
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PROPANE 


FIGURE  6.  THERMAL  NETWORK  FOR  A 
SINGLE  TUBE 
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illustrates  the  method  used.  If  UA  is  pictured  as  the 
reciprocal  ci  resistance  then  the  flow  cf  energy  fron  one 
fluid  tc  the  ether  can  be  pictured  as  a  series  electrical 
circuit  and  UA  is  determined  from 


OA  =  1/(R1  *  R2  +  R3  +  R4)  (21) 

Reference  15,  by  Holman,  shows  how  each  resistance  can  be 
found.  B1  is  the  thermal  resistance  due  tc  convecticn  on 
the  inside  cf  the  tubes,  the  sea  water  side. 


1 

H1  =  (28) 

7ThSW  dj  L 


R2  is  the  thermal  resistance  of  the  tube  wall. 


ln(d0/d,) 
H2  =  - ! —  (29) 

2  7T  kw  1 


R3   is   the   thermal  resistance  attributed  tc  boiling  en  the 
outside  cf  the  tubes. 


R3  =  (30) 


R4  is  the  thernal  resistance  due  to  any   other  resistances, 
such  as  fculing  and  corrosion  deposits. 


The  heat  transfer  coefficient  on  the  seawater  side, 
h.  ,  is  determined  from  the  dimensionless  group  called  the 
Nusselt  number  defined  as  follows  for  flow  inside  of  tubes. 


4-0 


Nu  =  -ill L  (31) 

Jc 


In  order  tc  iind  the  Nusselt  number,  three  distinct  regions 
of  flow  are  allowed  for  in  the  problem:  laainar, 
transitior,  and  turbulent  flow.  The  division  points  are  in 
terms  of  values  of  the  Reynolds  number,  Re  =  ud/3/H  * 

laninar:        0  <  Re  <  2,0C0 
transition:  2,000  <  Re  <  10,000 
turbulent:  10,000  <  Re  <  infinity 

In  the  laainar  region  the  Sieder-Tate  correlation  is  used. 


1/3 
NuL  =  1.86  (Re  Pr  dj /L  )  (32) 

In  the  turbulent  region  the   Eittus-Boelter  correlation   is 
used. 


0.8   0.3 

NuTu  =  0.023  Re   Pr  (33) 

In  the  transition  region  the  flew  is  unstable,  but  fron  data 

in  Hef.  14   it  appears  that  a  fairly  smooth  curve  drawn  from 

the   lamirar  curve   to  the   turbulent   curve  adaguately 

represents  the  value  of  the  Nusselt  number  in  the  transition 
region.  Ibis  is  done  by  assuming  the  form 


NuTr  =  y  NuL  +  (1  -  y)  NuTu  (34) 

where 


3      2 
y  =  afle  +  bRe  +  cRe  +  d  (35) 


41 


Th€  coefficients  in  the  polynomial  are  found  by  requiring 

NuTr  =  NuL  at  Re  «  2,000 
NuTr  =  NuTu  at  Re  =  10,000 

and   that  th€  first  derivatives  with  respect  to  the  Reynolds 
number  match  at  Re  equal  to  2,000  and  10,000. 

She  I.I.  Mostinski  correlation,  from  Ref.  16,  for 
single  horizontal  tube  pool  boiling  is  used  to  find  th€  heat 
transfer  coefficient  on  the  outside  of  the  tube.  BTBI's 
extensive  experimentation  led  Palen,  Yarden,  and  Tatcrek,  in 
Bef.  13,  to  conclude  that  Mostinski's  equation  produced  the 
"most  consistent  results".   His  correlation  is 


C.69  0.7       0.17      1.2       1C 
hwj  =  C.CC658  pc    q    (1.8  r     +  4  r    +  10  r   )  (36) 

where 

E  =  F/Pc 

pc  =  critical  pressure  of  the  working  fluid  (psia) 

g   =  heat  flux  (BTO/hr-f t2) 

A  final  heat  flux  relation  needed  is  the  relation  fcr  the 
maximum  heat  flux  allowed.  Mostinski,s  correlation  fcr  a 
single  tune  is  used. 

0.35       0.9 
9M«  =  803  Pc  r     <1  *  r>  <37> 

where  the  units  are  the  same  as  in  equation  (36) .  From 
HTEI*s  data,  under  unknown  conditions,  qmax  for  the  bundle 
aFfears  tc  te  one  third  of  the  qmQX  for  a  single  tube. 
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Ihe  power  required  to  pump  the  salt  water  ttrcugh 
the  tubes  depends  en  the  pressure  drop  across  the  tufce  bank. 
There  are  three  causes  of  the  pressure  drop.  The  flow 
experiences  a  loss  as  the  water  enters  into  and  exits  from 
the  tubes  in  addition  to  the  frictional  loss  along  the  tube. 
Frcm  Ref.  14,  the  pressure  drcp  fcr  a  liquid  is  determined 
frcm 


2    4  f  L 
Ap  =  0.5  p$w  uswv  (  _- ♦  Kc  -  Ke)  (38) 

di 


In  order  tc  determine  f,  Ke  and  Kc,  the  flow  regime  must  be 
specified.  The  regimes  are  a  little  different  than  those  for 
the  Kusselt  lumber. 

laminar:  0  <  Re  <  2,000 
transition:  2,000  <  Re  <  5,000 
turbulent:   5,000  <  Re  <  infinity 

The  Fanning  friction  factor,  f,  in  the  laminar  regime  is 

f  =  16/Re  (39) 

In  the  turbulent  regime,  the  Fanning  friction  factor  is 
found  frcm  the  Blasius  equation  in  the  range  5,O0C  <  Ee  < 
1CC,000 


0.25 
f  =  0.079/Re  (40) 


which   is   a   close   approximation  to  the  Karman  -  Nikuradse 
egcaticn, 
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-0.5  0.5 

<4  f)      =  -0.8  +  0.87  ln(Re(4  f)    )       <41) 

which  is  valid  throughout  the  turbulent  region  (Re  >  5,000), 
assuming  smooth  nailed  tubes.  Either  equation  (40)  or  (41) 
may  be  used  with  similar  results  but  the  Elasius  is  tetter 
suited  to  this  analysis  because  the  Karman-Nikuradse 
equation  must  be  solved  by  iteration.  Equation  (41) 
consumes  uore  computer  time  to  solve  and  can  cause  accuracy 
problems  if  not  enough  significant  digits  are  evaluated. 
The  fricticr  factor  in  the  transition  region  is  determined 
by  fitting  a  least  squares  curve  through  a  hand  drawn  curve 
from  the  end  of  the  laminar  curve  to  the  beginning  of  the 
turbulent  curve. 

The  entrance  effect,  Kc,  and  the  exit  effect,  Ke, 
depend  on  the  Reynolds  number,  the  free  flow  to  frontal  area 
ratio,  <t> ,  and,  in  the  laminar  region,  on  the  tube  length. 
For  the  rciler,  <f>  is  determined  by 


Q>  =   J (42) 

Ds   -  4  SA 


Figure  9,  Eage  7-13,  in  Ref.  17  is  used  to  find  Ke  and  Kc. 
The  effect  of  Reynolds  number  and  L/d  ratio  was  eliminated 
from  the  relationships  in  the  laminar  and  the  turbulent 
regimes  of  tte  flow,  tut  the  Reynolds  number  was  included  in 
the  traEsiticn  region  because  of  the  large  variation  of  Kc 
with  Reynolds  numter.  The  following  eguaticns  are  used  to 
find  the  entrance  effect  in  the  various  flow  regimes. 

Iaminai:  (using  the  curve  for  4  (L/d) Re  =  0.1) 

KcL  -  1.  -  0.4  4>  (43) 
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Turbulent:  (using  the  curve  for  Re  =  5,00C) 

KcTy  =  0.52  -  0.4  <j>  (44) 

Transition: 

KcTr  =  A  KcL  +  (1  -A)  KcTu  (45) 

where 

A  =  e  Re  +  f  (46) 

The  coefficients  e  and  f  are  determined  ty  requiring  the 
transition  Kc  to  egual  the  laminar  Kc  at  Re  equal  to  2,000 
and  the  turtulent  Kc  at  Re  egual  to  5,000.  The  exit  effect 
is  the  sane  fcr  all  flow  regimes,  and  is  found  frcn  the 
following  eouation  that  approximates  the  4  (L/d) Re  =0.05 
curve  in  Eef.  17. 


Ke  =  1.C0S3  -  2.5178<£  ♦  1.1613  <J>   -  0.17677(f)      (47) 


2-   ghe  Condenser 

The  design  of  the  condenser  is  almost  the  saae  as 
the  design  of  the  boiler.  It  is  a  fixed  tute  sheet, 
tube-in-sbell  heat  exchanger.  The  condenser  does  net  need  a 
separation  vclume  so  the  entire  shell  is  filled  with  tubes. 
In  order  to  allow  the  entering  vapor  to  reach  all  of  the 
tubes,   channels   through   the   tube   bank   are  needed.   The 
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allowance  for  this  is  accomplished  by  decreasing  the  shell 
diameter  by  a  fixed  factor  when  calculating  the  numter  of 
tubes.  The  calculation  of  the  working  fluid  side  heat 
transfer  coefficient,  h  ,  is  different.  The  equation 
accounts  for  the  condensate  from  one  tube  falling  en  those 
belcw  it  acd  for  the  subcooling  of  the  cendensate.  M.M. 
Chen's  eguaticn,  as  developed  in  Ref.  18  and  shown  in  a  more 
single  fcrm  in  Ref.  17,  is  used. 


3       „  0.25 


lwf 


where 


=    C.726 


h^Hh.a'1']    - 


k   =    BTD/ft2-hr-°F 

AT    =    TSflt     -    T„  (4  9) 

3 
H»       *   H,       +  —  cDAT  (50) 

n  =  number  of  tubes  in  a  column 

H«     is   the   enthalpy  corrected   for  subcooling   of   the 

f9 
ccrdensat€.   Eguaticn  (48)  replaces   eguaticn   (36)   ir  the 

analysis   cf   the   condenser.    In   the  circular  shaped  tube 

bark,  n  is  fcurd  by  determining  the  number   cf  tubes   in   a 

eclumn  of  average  height. 


Ds 
n  =  9.425  —  -  1  (51) 

Pi 


The  definiticc  of  AT-  becomes 


^6 


ATo  ■  **sot  -  Tsw>  (52) 

in  contrast  to  equation  (26)   above.    The   seawater   cutlet 
temperatuie  is  found  from 


-  T, 


AT, 


o 


which  replaces  equation  (25)  . 

In  cider  to  find  the  heat  transfer  coefficient,  the 
temperature  cf  the  outside  tube  surface  is  needed  in 
equation  (49)  ;  therefore,  the  thermal  network  is  solved  for 
the  outside  wall  temperature  and  for  the  heat  flow  rate. 

Q  =  UA»  (Tw  -  TbulK  )  (54) 

where 

UA»  =  1/(R1  ♦  R2  +  R4)  (55) 

Tw ii  =  average  of  the  inlet  and 
the  outlet  seawater  temp. 

There  is  nc  factor,  F,  in  the  heat  transfer  equation  for 
ccrdensation  since  the  effect  of  neighboring  tubes  is 
accounted  fcr  in  eguation  (47)  . 

3 .   11 e_ Eum£S 

Here  are  three  pumps  in  the  model.  Eoth 
circulation  pumps  for  the  boiler  and  the  condenser  are 
assumed   to  be  propeller  pumps  (axial  flow  pumps) .   The  feed 
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pumps  are  assumed  to  be  centrifugal  pumps.  The  differences 
in  the  characteristics  of  the  two  types  of  pump  show  up  only 
in  the  calculation  of  capital  costs  of  the  pumps  ir  the 
model. 

lie  chcice  of  the  type  of  pump  to  apply  is  made  on 
the  basis  cf  the  flow  rate,  GFM,  and  on  the  head  the  pump  is 
working  against.  Cualitatively,  the  propeller  pump  is  used 
where  the  flew  rates  are  high  and  the  head  is  low.  The 
centrifugal  pumps  apply  where  the  flow  rate  is  comparatively 
lew  and  the  head  pressure  is  high.  The  qualitative  picture 
is  placed  in  tetter  perspective  through  the  calculaticn  of 
the  specific  speed  from 


0.5 
BPM   GFM 
Nsg  =  (56) 

0.75 
Head 


The  centrifugal  pump  is  used  when  500  <  Nsg  <  7,500  and  the 
propeller  pu&p  applies  when  7,500  <  Nsg  <  15,000,  acccrding 
to  Bef.  1S.  The  choice  of  the  type  of  pump  is  checked  after 
the  optini2aticn  process  has  been  completed  by  solving 
eguation  (55)  .  If  the  value  of  the  specific  sp€ed  is 
outside  of  tie  range,  500  <  Nsg  <  15,000,  the  maximum 
efficiency  cf  the  pump  drops  off  rapidly;  therefore,  few 
pumps  are  designed  to  operate  outside  cf  this  range. 

Since  the  pumps  are  "black  boxes",  only  the  power 
consumed  by  the  pumps  is  calculated.  The  power  reguired  to 
dri^e  the  punp  is  found  from 


M  Ap 
P.  _L  ,57, 

fpumnp 
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E.   CCST  FEAHESOEK 


As  discussed  at  the  start  of  this  report,  the  capital 
costs  of  the  major  components  are  the  only  costs  to  be 
considered  it  this  mcdel.  The  boiler  and  the  condenser  are 
both  tute-in-shell  heat  exchangers  assumed  to  te  cf  the 
fixed  tube  sheet  type.  The  only  difference  being  that  the 
shell  of  the  boiler  is  partly  filled  with  tabes.  Eoth 
seawater  circulation  pumps  are  assumed  to  be  propeller  type 
puaps  and  nade  cf  the  same  material.  The  size  and  number 
per  heat  exchanger  shell  may  be  different  fcr  the  bciler  and 
the  condenser.  The  number  of  working  fluid  feed  pumps, 
which  are  the  centrifugal  type  pumps,  are  assumed  to  egual 
the  number  cf  bciler  shells. 

The  cost  estimating  relationships  fcr  the  heat 
exchangers  ard  the  centrifugal  pumps  are  taken  from  Bef.  20, 
by  K.fl.  Guthrie.  The  equation  fcr  the  cost  of  the  heat 
exchangers  is 


where  the  constants  are  calculated  from  a   leg-log   plct  of 
cost  versus  tctal  tube  surface  area  and 

Fj  =  0.6  for  a  fixed  tube  sheet 

Fh  =  0.C  for  pressure  correction  (p  <  150  psi) 

F^  =  material  factor  from  Table  1 

Ihc  =  ccst  index 

At  =    d0  L  N 

fcr  a  relevart  range  cf  100  <  A   <  10,000  ft*. 

The   ccst  eguaticn  for   the  centrifugal  pumps  and  the 
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associated  meters  is 


0.721 
*{p    =  1-C13  Fm  Fo  **p   CH  (59) 

where  the  cocstants  are  found  from  a  straight  line 
approximation  cf  the  log-log  plot  cf  cost  versus  C  for  the 
section  cf  tte  curve  from  C  =  30,000  to  C  ■  300,000  for  an 
electric  motcr  driven  pump  and 

F0  =  1.0  fcr  suction  temperatures  <  250°P  (121°C) 

I^p  =  ccst  index 

Fm  =  material  factor  from  Table  2 

Cm  ■  predect  of  the  flow  rate  in  GEM 

and  the  pressure  differential  in  psi 

for  a  relevant  range  of  30,000  <  CH  <  300,000.  The  cost 
data  in  Guthrie's  article  is  given  fcr  a  time  base  of 
oid-1968  . 

Reference  21,  also  by  Guthrie,  gives  cost  data  for 
propeller  funps  and  centrifugal  pumps.  The  eguation  fcr  the 
propeller  cuiips  with  the  motors  is 


0.783 
$cp    =  5.767  P^Fp  lcp    Flow  (60) 

where   the  ccrstants  are  found  from  zhe   log-leg  plot  cf  cost 
versus  flew  rate  for  an  electric  motor  driven  pump  and 


Fp  =  1.0  fcr  a  suction  pressure  <  150  psi 

F^  =  naterial  factor  from  Table  3 

ICp  =  cost  index 

Flow  =  flew  rate  in  GPM 

for   a   relevant   range  of  1,000  <  Flew  <  100,000  GEM.   The 
time  base  fcr  the  propeller  pump  data   is   fcr   the   end   of 
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Table  1 
MATERIAL  FACTORS  FOB  THE  HEAT  EXCHANGERS 


Material 

Shell 

Tube 

Fm 

Cartcn  Steel 

Carbcn  Steel 

1.00 

Careen  Steel 

Brass 

1.52 

Cartcn  Steel 

Stainless  Steel 

3.52 

Stainless  Steel 

Stainless  Steel 

4.50 

Cartcn  Steel 

Monel 

3.75 

Scnel 

Monel 

4.95 

Cartcn  Steel 

Titanium 

11.10 

Titaiium 

Titanium 

16.60 

Table  2 

MiTEflAL  FACTORS  FOR  THE  CENTRIFUGAL  FEED  PUEFS 


Material 

Fm 

Cast  Iron 

1.00 

Bronze 

1.28 

Cast  Steel 

1.32 

Stainless  Steel 

1.93 

Monel 

3.23 

Hastelloy  C 

2.89 

Titanium 

8.98 
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Table  3 
MATERIAL  FACTORS  FOR  PROPELLER  PUMPS 

pumps 

Material  Fm 

Cast  Iron  1 .00 

Cast  Steel  1.28 

Stainless  Steel  1.64 
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1S70.  Tie  costs  arrived  at  from  equations  (58),  (59),  and 
(60)  are  in  dcllars  per  component,  for  example,  dollars  per 
beat  exchanger  shell.  The  purpose  of  the  cost  index  is 
allcw  all  the  equipment  to  be  costed  on  a  ccmmon  basis.  The 
tine  difference  is  assumed  to  be  small  encugh  so  that  the 
cost  index  is  set  to  1.00  for  all  the  components. 


E.   ASSEKELI  CF  THE  MATHEMATICAL  MODEL 


The  engineering  and  the  ccst  relationships  are  linked 
together  in  the  model.  The  cost  equation  is  the  function  to 
be  minimized  and  the  implicit  and  explicit  constraints  are 
the  engineering  equations  that  restrict  the  feasible  region 
of  the  problem.  The  explicit  constraints  restrict  the 
variables,  the  vector  X,  to  positive  values.  The  rest  of 
the  model  is  fcmulated  in  the  following  form 

ninimize     f (X) 

subject  to:   g  (X)  >  0,  i=1,2,....,m 
i 

h  (X)  =  0,  i=m+1, . . . . ,m+n 
i 

where  the  definitions  of  the  elements  of  X  are  sho*n  in 
Table  4.  Several  of  the  variables  must  have  only  integer 
values  if  the  system  were  actually  built,  but  for  the 
objectives  cf  this  aralysis  they  are  left  as  continuous  real 
variables.  Icr  example,  the  number  of  shells  may  have  an 
integer  value  in  the  actual  system  but  this  requirement  can 
be  accomodated  by  solving  the  system  for  a  real  value  then 
integerizing  the  number  of  shells  and  solve  the  prcblen  with 
the  next  higher  and  the  next  lower  integer  number  of  shells. 
The  lower  cf  the  two  solutions  would  be  chosen.  The  final 
rescluticn   cf   the    correct    number    of    shells    would 
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Table  4 
DEEIKITICN  OF  THE  ELEMENTS  OE  THE  VECTOR  X 

x(1)  =  velocity  of  seawater  in  the  bciler 

i(2)  =    n  "     "     "    "   condenser 

i  (3)  =  cutlet  pressure  in  the  boiler 

x(4)  =  inlet     "     "   "   condenser 

x(5)  =  tube  length  in  the  boiler 

x(6)  *   ■     "    "    "  condenser 

x(7)  =   nass  flow  rate  of  the  working  fluid 

i(8)  =  diameter  of  the  boiler  shell 

x(9)  =    "     "    "   condenser  shell 

x(1C)  =  outside  dia.  of  a  tube  in  the  boiler 

2(11)  =    "     ■    "   "   "   "    "   condenser 

x(12)  =  number  of  bciler  shells 

x(13)  «   "    "   condenser  shells 
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prcbably  depend  on  other  than  thermodynamic  and  cost 
considerations,  such  as  a  requirement  that  one  extra  shell 
be  included  to  allow  for  one  shell  being  down  for 
maintenance.  Obviously,  the  number  of  tubes  in  a  shell  must 
be  an  irteger  value  but  one  extra  tube  cut  of  several 
thousand  is  a  minor  matter.  In  addition,  there  may  be  a 
requirement  fcr  extra  tubes  to  allow  for  a  percentage  to  be 
plugged  when  they  develop  leaks  so  that  the  entire  shell  is 
net  put  cut  cf  ccmnission  by  a  minor  leak  in  one  tube. 

1 •   £b jective  Function  and  the  Constraints 

The  objective  function  is  the  sum  of  the  costs  of 
the  boiler  ace  the  condenser,  and  the  pumps.  The  model 
assumes  that  there  is  one  circulating  pump  and  one  feed  pump 
for  each  shell  in  the  boiler  and  that  there  is  one 
circulating  pump  for  each  shell  in  the  condenser.  Since 
there  is  nc  cost  analysis  for  the  turbine,  nc  assumption  is 
needed  abcut  the  number  of  turbines  in  the  system. 

There  are  nineteen  constraints  in  the  present 
program,  cne  cf  which  is  optional.  The  first  three  cencern 
the  saturation  pressures  in  the  boiler  and  the  cendenser. 
The  next  twelve  place  restrictions  on  various  dimensions  of 
the  toiler  and  the  condenser.  The  last  three  require  that 
scue  intermediate  calculations  stay  positive. 

Csing  the  objective  function  as  a  constraint  has 
proved  advantacecus  in  cases  where  the  problen  was 
ncn-ccnvei  or  where  the  program  had  difficulty  in  finding  a 
solution.  It  is  used  by  setting  a  boundary  en  the  objective 
function  so  that  the  constraint  is  either  infeasible  or 
slightly  feasible  at  the  starting  point.  The  proper  choice 
must  be  found  by  trial  and  error.  One  caution  to  note  is 
when   equality  constraints   are   involved.   If  the  starting 
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point  is  near  the  optimum,  the  program  tries  to  move  away 
from  the  boundary  and  the  equality  constraints  may  not  be 
satisfied  when  the  program  finishes. 

The  need  for  the  ether  eighteen  inequality 
constraints  is  net  obvious  when  the  problem  is  first 
formulated.  The  thermodynamics  of  the  problem  require,  for 
example,  that  the  propane  saturation  temperature  be  less 
than  seafcater  temperature  in  crder  for  the  toiler  to  absorb 
energy,  tut  the  program  has  no  way  of  knowing  that,  unless 
it  it  specified.  If  the  program  can  improve  the  cost  and 
satisfy  the  ether  constraints  by  reversing  the  relationship 
it  will  dc  so.  This  inversion  may  take  place  only  for  a 
short  tine  during  the  solution  of  the  problem,  but,  if  the 
equations  fail  under  this  condition,  the  constraint  en  the 
temperatures  in  the  boiler  is  necessary.  Most  of  the  other 
constraints  are  needed  for  the  same  reason. 

Seme  constraints  are  not  added  to  the  problem  since 
they  are  net  binding  in  the  solution.  There  are  twe  such 
censtraints.  The  first  is  the  segmental  area  required  to 
allow  separation  of  the  liquid  droplets  from  the  vapor.  The 
segmental  area  set  aside  by  the  40£  rule-cf-thumb  is  more 
that  that  required  by  the  consideration  of  the  physical 
properties  cf  prcpane.  The  ether  constraint  is  the 
restriction  en  the  maximum  heat  flux  by  Mostinsk^s 
correlaticn,  equation  (37).  The  constraints  dc  not  become 
binding  in  the  solution. 

The  three  equality  censtraints  contain  the 
calculations  cf  the  properties  of  the  various  components 
modeled  in  this  analysis.  The  first  cne  concerns  the 
boiler.  It  states  that  the  actual  energy  absorption  rate  of 
the  boiler  shells  Bust  be  equal  to  the  energy  absorption 
rate  reguired  by  the  thermodynamic  cycle.  Similarly,  the 
second  eguality  constraint  states  that  the  heat  rejected   by 
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the  condenser  must  equal  that  required  by  the  thermodynamic 
cycle  for  the  condensation  of  the  vapor  from  state  point  5 
back  to  state  pcint  1  as  shown  in  Fig.  3.  The  third 
constraint  calculates  the  power  required  by  the  various 
pumps  and  the  pcwer  extracted  from  the  turbine  to  run  the 
system*s  pumps  and  to  produce  electrical  power  for 
consumption.  This  constraint  requires  that  the  power  cut  of 
the  turbine  equal  the  power  used  by  the  pumps  and  the 
generator.  The  net  useful  power  cut  of  the  systen  is  a 
design  objective  and  is  specified  as  25MW  herein. 

2-   The  Optimization  Method 

Since  the  arrival  of  the  electronic  computer,  a 
great  many  nonlinear  programming  technigues  have  been 
developed  frcm  early  theories  and  more  current  research. 
The  computer  program,  called  the  Sequential  Unconstrained 
Minimization  Technigue  (SOMT) ,  developed  by  Mylander, 
Helens,  and  Mccormick  in  fief.  22  is  used  to  optimize  the 
model.  Their  program  implemented  much  of  the  theory 
contained  in  fief.  23,  by  Fiacco  and  McCcrmick,  about 
ncnlinear  programming  using  unconstrained  minimization 
technigues.  SDMT-Version  4  is  chosen  because  of  its  ability 
to  solve  a  wide  variety  of  problems.  The  theoretical 
requirements  for  SOMT  to  be  guaranteed  to  find  a  local 
minimum  is  that  there  exists  seme  point  which  satisfies  the 
ineguality  constraints  and  that  there  not  be  a  local  minimum 
at  pcints  where  X  goes  to  infinity.  If  certain  convexity 
conditions  are  satisfied,  then  the  local  irinimum  found  is  a 
global  minimum. 

SCMT  uses  a  penalty  function  to  transform  the 
constrained  optimization  problem 
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ninimize     f (X) 

subject  tc:   g(X)  >  0,  i=1,2, ,m 


1 


h  (X)  =  0,  i=m+1 , . . . . ,m+n 
i 


into  a  seguence  of  unconstrained  problems  cf  the  form 


P(X,r)  =  f  (X)  -  r.I  In  g  (X)  +  1/r  I   (h  (X) )  2    (61) 

i=1    1         i=m+1   i 

where  E(X,r)  becomes  the  function  to  be  minimized.  The 
inequality  constraints,  g  (X)  ,  cause  a  large  penalty  to  occur 
as  a  boundary  is  approached.  The  equality  constraints  add  a 
large  penalty  whenever  h (X)  departs  from  zero  in  either 
direction.  The  size  of  the  penalties  due  tc  g (X)  and  h (X) 
is  controlled  by  the  arbitrary  parameter  Hr".  The  algorithm 
solves  the  problem  by  solving  a  series  of  surproblems.  For 
r1 >r^>rg . . .>0,  SOMT  minimizes  P(X(rrt),rn)  in  each  sutproblem 
in  order  to  find  the  minimum  of  f  (X)  as  "r"  tends  tc  zero. 
When  a  minimum  of  a  subproblem  is  found,  "rM  is  reduced  by  a 
factor  assigned  by  the  user  and  the  new  subproblem  is 
solved.  The  process  continues  until  the  convergence 
criteria  chosen  by  the  user  is  satisfied.  Figure  7 
illustrates  schematically  the  results  cf  the  sequence  of 
subprcbleis  and  the  qualitative  shape  cf  the  penalty 
function.  Three  methods  of  minimization  are  included  in  the 
program.  The  user  may  choose  either  the  generalized 
Newton-Rapbscn,  steepest  descent,  or  McCormick  modification 
of  the  Fletcher-Powell  method. 

There  are  several  other  options  that  control  the 
operation  of  SOMT.  The  user  may  cheese  the  initial 
penalty,  1,  cr,  if  the  problem  satisfies  certain  conditions, 
the  user  may  elect  to  have  SOMT  compute  the  initial  "r". 
There  are  three  choices  for  the  convergence  criterion  cf  the 
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FIGURE  7.  PROBLEM  SOLUTION  SEQUENCE 
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prcblem.  The  criterion  used  depends  on  the  ineguality 
penalty  tera  teing  close  enough  to  zero.  Another  option  is 
the  choice  of  the  calculation  of  the  sutproblem  stepping 
criterion  .  Two  of  the  choices  compute  the  inverse  cf  the 
secend  partial  matrix  and  the  third  choice  computes  only  the 
gradient  of  the  penalty  function  (VJfP(I,r)).  One  cf  the 
mest  useful  cptions  is  extrapolation.  After  two  or  three 
sutproblems  have  been  solved  SUMT  is  able  to  make  a  first  or 
secend  crcer  extrapolation  to  the  starting  pcint  for  the 
next  sufccrcfclem.  This  can  speed  up  the  solution  cf  the 
problem  measureably. 

The  characteristics  of  the  problem  limit  the  choices 
for  each  cpticn.  The  penalty  factor  couldn*t  be  set  by  SDMT 
due  to  the  presence  of  equality  constraints  and  because  the 
initial  1  is  usually  not  close  to  any  boundary.  The  stepping 
criterion  fcr  the  problem  is  that  the  ineguality  penalty 
ten  beccme  less  than  some  value,  usually  10-7. 
Extrapolation  is  used,  but  it  makes  little  difference  in  the 
answer  or  to  the  solution  time  whether  a  twe  point  or  three 
ccint  estimate  is  made.  The  subproblem  convergence 
criterion  is  the  magnitude  of  the  gradient  of  E(X,r) 
beccming  less  than  some  small  number. 

The  choice  cf  starting  penalty,  r,  affects 
significantly  the  solution  time  and  even  whether  a  solution 
is  found  ir  seme  problems  according  to  Bef.  7.  If  the 
initial  "r"  is  toe  small,  the  subproblem  optimum  may  te  the 
prcblem1 s  cptimum  which  may  be  hard  to  find  cr  require  more 
calculaticns  than  if  a  much  larger  "r"  is  chesen.  If  "r"  is 
toe  large,  the  first  few  optima  lie  near  the  center  of  the 
feasible  region  and  are  unrelated  to  the  optimum  cf  the 
objective  function. 

Scaling  of  the  variables,  constraints,  and  the 
objective  function  is   necessary.    The   variables   must   be 
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within  two  cr  three  orders  of  magnitude  of  each  other.  For 
example,  the  working  fluid  flow  rate  is  of  the  order  of  107 
whereas  the  tube  diameters  are  of  the  order  10°.  If  the 
finite  difference  interval  is  10-*,  the  partial  of  the 
eguations  with  respect  to  the  mass  flow  rate  may  be  lost  due 
tc  the  roundoff  error  providing  the  only  significant 
figures.  Hith  the  constraints  and  the  objective  function, 
unsealing  nay  result  in  one  term  in  the  penalty  furction 
being  dcninant,  and  the  algorithm  finds  that  most  cf  the 
decrease  in  the  penalty  function  takes  place  in  the  dcainate 
term.  The  result  is  that  SOMT  works  hard  to  decrease  that 
term  and  does  not  see  the  other  terms  resulting  in  a 
meaningless  sclution. 

In  making  first  and  second  order  estimates  of  the 
scluticn,  SCAT  may  generate  negative  numbers  in  the  X 
vector.  Hhen  it  evaluates  the  constraints  and  objective 
function  for  this  point,  the  eguations  fail  and  may  cause 
the  program  tc  fail  depending  on  the  severity  of  the  failure 
and  how  many  times  the  error  has  cccured.  This  problem  is 
circumvented  ty  testing  the  vector  upon  entry  intc  the 
subroutine  where  the  objective  function  and  the  constraints 
are  evaluated  and  setting  the  constraint  being  evaluated 
egual  to  seme  negative  constant. 

The  main  program  of  SUMT  is  altered  so  that  the 
prcblem  can  be  restarted  at  the  current  solution  without 
reading  in  a  new  problem.  This  is  necessary  because  the 
prcgram  would  give  all  the  indications  of  having  fcund  a 
sclution,  tut  if  the  problem  is  restared  it  could  oftec  find 
a  lew  value  for  the  cbjective  function.  The  new  soluticn  is 
compared  to  the  previous  solution  and  if  the  two  solutions 
are  net  clcse  enough  the  problem  is  restarted  after  aaking 
seme  changes  to  particular  parameters,  options,  or 
constraints  depending  on  the  objective  the  user  has  in  mind. 
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III.   RESULTS  AND  CONCLUSIONS 


One  cf  the  majcr  tasks  in  th€  development  of  a  nonlinear 
model  is  the  linking  of  the  model  with  the  optimization 
program  in  order  to  analyze  the  model.  The  selection  cf  the 
proper  nonlinear  programming  technique  is  essential  to  the 
solution.  Several  techniques  require  that  the  objective 
function  and  the  constraints  be  in  a  particular  form.  For 
example.  Geometric  Programming  developed  by  Euffin,  Peterson 
and  Zener,  in  Ref.  24,  requires  that  the  equations  be 
posyncmials  which  is  not  true  of  the  model  developed  in  this 
paper.  SOMI  was  chosen  because  cf  its  wide  applications  to 
problems  cf  any  type.  It  only  reguires  that  there  exist 
some  feasible  point  and  that  local  minima  not  exist  at 
infinity.  It  has  successfully  solved  many  ncnconvex 
problems. 

Two  versions  of  the  analytical  model  are  presented.  The 
first  is  called  the  iterative  version  because  the 
calculation  cf  the  heat  absorption  rate  and  heat  rejection 
rate  reguire  iterating  on  some  quantity  to  find  them. 
Mostinski's  correlation,  equation  (36) ,  contains  the  heat 
flux  which  is  an  unknown.  It  is  not  possible  to  solve  the 
system  cf  equations  explicitly  for  the  heat  flux;  therefore, 
a  successive  approximation  technique  is  used.  The 
calculation  cf  the  heat  rejected  in  the  condenser  involves 
the  same  process  with  the  temperature  of  the  outside  tube 
wall  as  the  iterating  variable.  SDMT  is  urable  to  arrive  at 
a  consistent  solution  in  the  time  allowed,  up  to  90  minutes 
of  computer  time.  The  solutions  are  quite  different  for 
each  try  vith  nc  pattern  evident  about  the  direction  the 
solution  might  lie. 
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There  are  several  possible  explanations  fcr  this 
behavior.  Perhaps  the  behavior  of  this  version  cf  the 
problem  is  such  that  SDMT  is  not  the  best  choice  fcr  the 
optimizing  algorithm  or  possibly  the  prctlem  cannct  be 
sclved  at  all  in  this  form.  Another  source  of  prcblems 
could  be  the  numerical  differentiation  procedure.  Numerical 
differentiation  is  subject  to  severe  truncation  and 
rcund-off  error.  The  subproblem  always  terminates  when  SUMT 
is  unabl€  to  reduce  the  penalty  function  during  the  next 
iteration,  frcm  one  estimate  of  X  to  the  next  estimate. 
Remember  the  convergence  criterion,  never  satisfied  for  this 
version  cf  the  problem,  is  the  magnitude  cf  V^E(X,r) 
becoming  less  than  some  small  number.  No  conclusive  results 
could  te  obtained  with  this  version  so  the  second  version  is 
used  fcr  the  present  analysis. 

In  this  second  version  the  heat  transfer  coefficierts  on 
the  Eropane  side  of  the  tube  are  fixed.  This  enables  the 
eguality  constraints  for  the  boiler  and  the  condenser  to  be 
calculated  to  the  full  accuracy  of  the  computer  without 
increasing  the  time  reguired.  The  correlations,  eguation 
(36)  and  (48) ,  are  used  for  comparison  with  the  fixed  ?alue. 
If  the  calculated  value  and  the  fixed  value  are  within  20  to 
2531  of  one  another,  the  solution  would  be  considered  to  be 
gccd  sirce  this  is  the  tolerance  claimed  for  the 
ccrrelaticns. 

SDMT  seemed  to  be  able  to  solve  the  new  prcblera  because 
a  series  cf  seven  solutions  have  been  produced  where  the 
values  of  the  elements  of  X  and  cf  the  objective  furction 
vary  nc  more  than  3.7%  and  1.2*,  respectively,  about  the 
mean  when  different  starting  points  and  different  SUMT 
control  parameters  are  selected.  See  Table  6  for  a  summary 
of  the  solutions  ottained  from  the  assumed  data,  as  shewn  in 
Table  5.  The  design  implications  of  these  results  are  guite 
reasonable.   The  range  shows  the  spread   in   the   solutions. 
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Table  5 
PROBLEM  INPUT  DATA 


Fouling  Thermal  Besistance 
Eciler  and  Condenser 

Elevation  cf  Boiler 
Afccue  Condenser 

Channel  Allowance 

Efficiencies  for  the  Turbine 

and  the  Pumps 
Seawater  Temperatures  -  Hot 

-  Cold 

Set  Output  Pcwer 

Eundle  Jactor  (F) 
Cost  Icdexes  (I) 
Tube  Rail  Thickness 

Min.  Distance  Between  Tubes 

Heat  Transfer  Coefficients 
fcr  Eciler  and  Condenser 


0.0  hr-0F/BTO 
(0.0  °C/kW) 
15  ft 
(4.57  m) 
0.83  ft 
(0.25  m) 
8551 

75°F 
(23.9°C) 
40°F 
(4.4<>C) 

8.532x107  BTU/hr 
(25  MS) 
3.0 

1 .00  fcr  all  ccmp 
0.028  in 
(0.071  cm) 
0.5  in 
(1.27  cm) 

145  BTG/ft2-hr-°F 
(0.822  kW/m2-oC) 


Materials 


Eoiler  and  Condenser 

Feed  Pumps 
Circulation  Jumps 


Carbon  Steel  Shells 
Brass  Tubes 
Bronze 
Stainless  Steel 
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Table  6 
EXAHPLE  P50BLEU  SOLUTION  SUHHAEI 


Variable 

Range 

Mean 

Var 

~2L 

Units 

Bciler : 

Velocity 

3.323-3.347 

3.33  5 

0.7 

ft/sec 

(1.013-1.020) 

(1.016) 

(a/sec) 

Pressure 

ajprox  0 

118.2 
(815.0) 

0.0 

lb/in* 
(*P) 

Tuhe  length 

39.16-40  .00 

39.71 

2.1 

ft 

(11.91-12.19) 

(12.10) 

(m) 

Shell  dia. 

24.12-25.00 

24.91 

3.5 

ft 

(7.352-7.620) 

(7.593) 

(m) 

lube  dia. 

0.756-0.761 

0.758 

0.7 

in 

(1.919-1.934) 

(1.926) 

(cm) 

#  of  shells 

7.  141-7.409 

7.228 

3.7 

- 

Ccrdenser : 

Velocity 

2.793-2.823 

2.808 

1.1 

ft/sec 

(0.851-0.860) 

(0.856) 

(ra/sec) 

Pressure 

91.53-91.60 

91.56 

0.1 

lb/in2 

(631.  1-631.6) 

(631.3) 

(*P) 

lube  length 

39.35-40  .00 

39.78 

1.6 

ft 

(11.99-12.19) 

(12.12) 

(m) 

Shell  dia. 

44.65-45.00 

44.89 

0.7 

ft 

(13.62-13.72) 

(13.68) 

(m) 

lube  dia. 

0.7*32-0.7686 

0.7660 

0.7 

in 

(1.938-1.952) 

(1.946) 

(cm) 

I  of  shells 

3.297-3.401 

3.331 

3.  1 

- 

Mass  flc«  rate 

22.12-22.21 

22.17 

0.4 

10°lbra/hr 

(10.03-10.07) 

(10.06) 

(10*kg/hr) 

System  ccst 

7.578-7.673 

7.607 

1  .2 

106  3 
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The  arithmetic  mean  of  the  seven  solutions  is  calculated. 
The  difference  between  the  maximum  and  the  minimum  values  is 
shewn  as  a  percent  of  the  mean.  The  percent  variation  gives 
seme  idea  cf  how  well  the  solution  was  determined.  Note 
that  the  tufce  length  and  shell  diameter  constraints  are 
binding  for  both  the  boiler  and  the  cendenser.  The 
sutFrcblem  still  terminates  only  when  the  penalty  furction 
does  not  decrease  from  one  iteration  to  the  next.  The 
magnitude  cf  \^P(X,r)  is  not  becoming  small  enough  to 
satisfy  tte  convergence  criterion.  The  value  never  becomes 
less  than  C.02  and  is  usually  greater  than  0.5  *hen  a 
sutprcbleir  is  stopped  by  the  program. 

Useful  results  are  obtainable,  however.  The  engineer 
can  learc  much  about  his  problem  even  from  a  single  point. 
Fcr  exaaple,  from  the  cost  breakdown  iD  Figure  €r  the 
portion  cf  the  costs  attributable  to  each  component  can  be 
seen.  Ihe  heat  exchangers,  of  course,  dominate  with  785?  of 
the  cost.  Tte  circulation  pumps  are  about  219?  of  the  cost 
with  the  feed  pumps  being  T.2J5  of  the  total  cost.  Since 
the  efficiency  assumed  for  the  circulation  pumps  is  prcbafcly 
toe  high,  these  pumps  would  be  even  mere  ccstly  in  terms  of 
energy  and  meney.  This  could  indicate  the  priorities  that 
should  re  set  on  the  design  efforts  for  improvements  cf  the 
vaiicus  components  with  the  highest  priority  being  the 
highest  percentage  cost  component. 

The  systems  pewer  distribution  is  shown  in  Table  7. 
The  significant  features  are  the  large  energy  absorption 
rate  for  the  snail  useful  output  power  and  the  parasitic 
pumping  pewer.  The  lower  the  efficiency  the  greater  must  be 
the  power  absorbed  for  the  same  output.  This  is  to  be 
expected  fcr  this  kind  of  system.  The  parasitic  pumping 
pewer  (12*  of  the  turbine  output)  is  about  twice  what  it  is 
for  a  fessile  fuel  plant. 


66 


t   CONDENSER 

V- CIRCULATION 

PUMPS 

FEED     / 

BOILER 

$978,000 

PUMPS    t— 

—  CIRCULATION 

12.9% 

$94,100 

PUMPS 

1.2* 

$654,000 
8.61 

FIGURE  8. 

COST  BREAKDOWN 
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Table  7 
SOLUTION  POWEE  EUEGET 


Eeat  Absorption 

Turbine  Shaft  Power 

Jeed  Pumps 

Elr  Circulation  Pumps 

Ccnd  Circulation  Pumps 

E€at  Eejected 

Ket  Cseful  Power 


BTO/br 

MWatts 

3.56x1C9 

1042.8 

9.69x10^ 

28.4 

4.49x106 

1.3 

2.58x10* 

0.8 

4.51x10* 

1  .3 

3.47x1C* 

1015.7 

8.53x107 

25.0 
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Once  the  model  has  been  optimized  with  the  number  of 
shells  as  real  variables,  the  problem  can  te  re-solved  twice 
by  fixing  the  shells  at  the  next  higher  and  the  next  lower 
integer  cutter.  For  the  higher  number  of  shells,  the  cost 
increased  slightly  as  expected,  from  7.607  million  dcllars 
average  to  7.983  million  dcllars.  The  problem  with  the 
lower  integer  number  of  shells  could  not  be  solved  due  to 
the  constraints  en  the  tube  length  and  the  diameter  cf  the 
shells  feeing  binding  before  the  shell  numbers  were 
integerized.  The  program  is  unable  to  find  a  new  solution 
that  satisfies  the  eguality  constraints. 

Table  8  shows  some  miscellaneous  data  from  the  model 
that  deserves  comment.  The  overall  heat  transfer 
coefficient  cf  the  boiler  is  approximately  three  times  that 
for  the  cendenser  which  is  the  major  reason  for  the  lewer 
ccst  of  the  boiler  relative  to  the  condenser.  However, 
rememher  that  the  bundle  factor,  egual  to  3,  is  a  part  of 
the  heat  transfer  coefficient.  More  data  must  be  gotten 
abcut  this  factor  in  crder  to  improve  the  reliance  oe  any 
solution.  There  is  a  significant  difference  between  the 
Prcpane  heat  transfer  coefficients  that  were  assumed  and  the 
theoretical  values  calculated  from  the  correlations.  Other 
values  shculd  be  tried  until  the  fixed  value  is  withir  the 
209  tolerance  on  the  theoretical  value.  The  flew  of 
seawater  through  the  tubes  in  both  the  boiler  and  the 
ccrdenser  is  turbulent.  Apparently,  the  resulting  cost 
savings  due  to  the  improved  heat  transfer  rate  through  a 
smaller  area  is  more  than  the  increased  costs  of  the  higher 
pumping  power  needed.  It  seems  likely  that  the  solution 
obtained  is  globally  optimal;  several  different  starting 
pcints  were  tried.  The  problem  has  been  started  in  both  the 
laiinar  and  the  turbulent  flow  regions  with  no  change  in  the 
results. 


69 


Table  8 
ADDITIONAL  DATA  FROM  SAMELi 
PROBLEM  SOLUTION 

Boiler   Condenser  Oni*s 

Scrking  fluid  Heat  Transfer  Coefficients 

assumed                  145.        145.  BTO/hr-f t*-°P 

(0.822)     (0.822)  (kW/n2-cC) 
Theoretical             227.       241. 

(1.29)      (1.37) 

Surface  ir«a  per  Shell   2.39x10*    1.27x10*  ft2 

(2.22x10*)  (1.18x105)  (m2) 

Seawater  Cutlet  Temp       72.3       43.8  °F 

(22.4)      (6.55)  (°C) 

Shell  Saturation  Temp      66.7       49.9  °F 

(19.3)      (9.94)  (°C) 

Seawater  Bass  Plow  Rate  4.39x10s    9.11x1fl«  lbm/hr 

(1.99x10*)  (4.13x10^)  (kg/hr) 

Mean  Temp  Difference       6.67       7.85  °F 

(3.71)      (4.36)  (°C) 

Overall  Beat  transfer     310.6      105.0  ETO/hr-f t2-°F 

Coefficiert            (1.76)     (0.595)  (kW/ffl2-oC) 

Heat  Flui                2070.0       824.9  ETU/ft2-br 

(6.528)     (2.602)  (kW/m2) 
Seawater  Beyrclds  eo.    19,800.     9,960. 
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IV.   RECOMMENDATIONS 

The  present  model  has  a  great  deal  mere  inforaati.cn  to 
offer,  but  seme  improvements  in  the  model  should  be  made  to 
improve  the  ability  of  a  nonlinear  program  to  solve  the 
prcblem.  The  first  recommendation  would  be  to  stud?  the 
equality  constraints  and  place  inequality  constraints  en  any 
variable  cr  grouping  of  variables  if  they  fceth  can  tecome 
negative  and  also  must  have  their  logarithms  computed. 
Examples,  already  included  in  the  program,  are  the  numter  of 
tures  in  the  shells.  The  heat  transfer  coefficients  fcr  the 
working  fluid  side  in  the  boiler  and  the  condenser  could  be 
included  as  equality  constraints.  This  can  be  dene  by 
adding  the  heat  transfer  coefficients  to  the  vector  X  and 
requiring  that  the  variables  be  equal  to  the  coefficierts  as 
cemputed  ry  the  appropriate  correlation.  The  same  could  be 
dene  with  Ealen  and  Smalls  property-dependent  calculation 
of  the  separation  area  by  equation  (12)  and  (13).  Frcm  an 
overall  standpoint,  the  results  seem  to  indicate  either  the 
prcblem  is  tec  difficult  to  be  solved  efficiently  ir  its 
present  fori  or  that  another  nonlinear  prcgram,  other  than 
SOMI,  might  te  better  suited  to  the  problem.  Numerical 
differentiation  of  the  constraints  in  this  problem  do  not 
seem  to  help  the  program  to  find  the  solution  since  the 
magnitude  of  \  P  (X,r)  is  not  near  zero  when  a  sutprotlem 
steps. 

Once  the  optimization  scheme  is  functioning  properly, 
the  model  can  te  used  to  study  the  effects  cf  fouling  and 
ccrrosicr  and  the  costs  of  their  prevention.  Many  proposals 
have  teen  made  that  must  be  analyzed  in  the  ccntext  of  the 
total  sjsten.  Seme  such  as  Amertap,  where  plastic  fcam 
tails  are   pumped   through   the   tubes   to   clean   the   heat 
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transfer  surfaces,  add  a  cost  penalty  from  the  capital  and 
operating  cost  of  the  Amertap  system  and  from  the  increased 
pumping  power  required.  Coatings  add  the  capital  costs  of 
the  coatirg  and  the  added  capital  cost  of  increased  heat 
transfer  area  reguired  because  of  the  increased  thermal 
resistance  cf  the  coating. 

A  sensitivity  analysis  can  be  performed  on  the  model  to 
test  the  assumptions  and  the  various  constants.  The 
sensitivity  to  changes  in  these  factors  can  indicate  how 
much  one  should  be  willing  to  pay  for  the  improvement 
sought.  High  sensitivity  indicates  areas  where  the  greatest 
effort  for  improvement  should  be  desired  or  the  priorities 
shculd  be  set  en  research  into  the  unknown  areas  cf  the 
problem.  The  trade-off  between  cost  and  performance  can  be 
studied  by  changing  the  output  power  and  looking  fcr  the 
lowest  ccst  per  kilowatt-hour,  which  is  the  measure  by  which 
ccnventicnal  pewer  plants  are  compared.  One  could  lcok  into 
the  cost  and  the  benefits  of  increasing  the  depth  from  which 
the  cold  water  is  pumped  versus  the  improved  temperature 
differential  possible  from  using  the  colder  water. 

For  the  capital  cost  of  the  system  to  be  complete,  the 
turbine  design  must  be  included  and  ccst  estimates  made. 
The  turbine  will  be  unique  due  to  the  very  low  pressure 
potential  drcp  it  must  use.  Some  work  has  been  done  by 
ether  grcups  working  on  OTEC,  Bef.  25  and  26.  Reference  27, 
by  Bcbert  L.  Bartlett,  can  provide  costing  and  design 
infermatien  although  he  deals  only  with  steam  turbines.  The 
remaining  part  of  the  system  that  is  desireable  for 
inclusion  in  the  analysis  presented  herein  is  the  cold 
seawater  intake  pipe.  It  will  be  on  the  order  of  1,0C0  ft 
(3C4.6m)  leng  and  100  ft  (30.4  m)  in  diameter. 
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